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Abstract 

Based on the formahsm developed by Psarobas et al [Phys. Rev. B 62, 278(2000)], which using 
the muhiple scattering theory to calculate properties of simple phononic crystals, we propose a very 
simple method to study the NaCl-type phononic crystal. The NaCl-type phononic crystal consists 
of two kinds of non-overlapping elastic spheres with different mass densities, Lame coefficients and 
radius following the same periodicity of the ions in the real NaCl crystal. We focus on the (001) 
surface, and view the crystal as a sequence of planes of spheres, each plane of spheres has identical 
2D periodicity. We obtained the complex band structure of the infinite crystal associated with 
this plane, and also calculated the transmission, refiection and absorption coefficients for an elastic 
wave (longitudinal or transverse) incident, at any angle, on a slab of the crystal of finite thickness. 
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I. INTRODUCTION 



Recently, inspired by the remarkable properties of the photonic crystals^, there has been 
increasing interest in the study of phononic crystals^'^'^'^i^, which are constructed by compos- 
ing identical particles (the inclusions) on host medium periodically. The acoustic behaviors 
in such materials are different from those in ordinary materials, such as the frequency band 
structure of classical waves and the propagation of elastic waves in the materials. Various 
applications are also proposed, for example, acoustic lensii, waveguides^ and negative re- 
fraction materialgi^. The behaviors of the new materials are related to their band structures 
so that the calculation of the acoustic wave spectrum, or band structure is the key point in 
the studies of phononic crystals. 

There are several techniques for the calculation of frequency band structures, the plane- 
wave approach is the simplest though the convergence is slow^i^*^, the multiple-scattering 
theory (MST) or KKR (Korringa, Kohn, and Rostoker) approach are widely used with 
better convergence^'^'i^, and the finite difference time domain (FDTD) method is a new and 
promising method^^. When the inclusions of phononic crystals are spheres in 3-dimensional 
case or cylinders in 2-dimensional case, layer MST is a very effective calculation technique 
not only for the evaluation of band structures, but also for the calculation of the transmission, 
reflection and absorption coefficients of a slab of the crystal of finite thickness. Based on 
the layer MST approach, Psarobas et al^*^ have developed a formalism for calculating the 
frequency band structure of phononic crystals and the corresponding scattering properties. 
In their algorithm and the related program package, they consider only the 3-dimensional 
phononic crystals formed by one inclusion per unit cell, so the phononic crystals are the 
so called "simple" crystals. In addition to band structure, their program package can also 
deal with scattering problems with a slab of crystals of finite thickness, if all successional 
scattering planes are 2-dimensional "simple" lattices. 

However, it is desirable to have a technique to deal with "complex" phononic crystals (i.e. 
crystals have more than one type of inclusion per unit cell). The crystals with "complex" 
structures may offer us different band gaps and new properties which are not existed in 
"simple" crystals. By generalizing the simple lattice layer MST method, we propose here 
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a simple technique to deal with NaCl-type "complex" phononic crystals by means of MST. 
Calculation of the band structures is made and compared to the "simple" case. A general 
method to work with "complex" phononic crystals, which are similar to those in low energy 
electron diffraction of "complex" crystals^^, is also discussed in Appendix A. Our method is 
based on MST and different from the plane- wave approach suggested by Zhang et al^^, so it 
is more efficient in dealing with spherical inclusion, and the case of fluid-solid composites. 
In addition, one may notice that the NaCl-type crystal can also be seen as a succession of 
two nonoverlapping planes of different spheres in the (111) direction and considered as a 
heterostructure^ or as a case of a super cell with planar defects^, so the complex band struc- 
ture could be extracted from their program. Actually, Liu et al^ have also used their bulk 
KKR method to calculate the full band structure of some complex phononic crystals with 
HCP and Diamond structure. However, the layer KKR method here we used is necessary, 
since it can also be used to calculate the transmission, reflection and absorption coefficients 
for an elastic wave incident on a slab of NaCl-type "complex" crystal. Recently, it has also 
been applied to a chiral structure, which is regarded as a more complex heterostructure, 
many interesting results was found^'^. The programm package developed by Sainidou et aU^ 
is a stable and efficient package in treating the simple phononic crystals, we have used the 
package extensively in our "complex" phononic crystal calculation. 

The paper is organized as follows. In section [Tll we present a simple method to evalu- 
ate the scattering properties of a slab consisting a number of layers, which could contain 
two kinds of spheres with specific 2D periodic structure. In section IIIH we demonstrate 
the applicability of our method with some numerical results and discussions. Finally, the 
conclusions of the paper are given in Sec. [IVl 

II. THEORY 

In this section, we briefly review the theory of layer MST used in the calculation of 
band structures of phononic crystals with the NaCl-type structure. Firstly, we consider a 
"complex" 2D crystal in xy-pla.ne at z = 0. There are two kinds of spheres in the 2D crystal. 
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The first is denoted as type A, located on the sites of a 2D lattice specified by 

R„, = niHi + ?7,2a2, (1) 

where ai,a2 are the primitive vectors of a 2D lattice in the xy-plane, and ni,n2 = 
0, ±1, ±2, ±3, .... The second kind of spheres is denoted as type B, located at 

= ab + R„ = ab + riiai + n2a2, (2) 

where ab = (ai + a2)/2 for NaCl-type phononic crystal. In this case, there are two types of 
inclusion per unit cell, and the corresponding 2D reciprocal lattice is 

g = mibi + m2b2, (3) 

where mi, m2 = 0, ±1, ±2, ±3, ... and vectors bi, b2 are defined by 

hi-aj = 27rSij, z,j = l,2. (4) 

Now we consider that a plane wave, either longitudinal or transverse, is incident onto 
such a plane of spheres of type A and B. The displacement vector of the incident wave has 
the forn>^, 

Uil(r) = "^[uinji'i' exp(iKj^, ■ r)ei,, (5) 

i' 

where s' = +(— ) represents that a wave is incident onto the plane from the left(right), u' 
denotes the polarization of the incident wave: q^' = qi = uj/ci for a longitudinal wave and 
qu' = qt = ^/ct for a transverse one. Because of the 2D periodicity in the plane of spheres, 
wavevector can be rewritten as 

Kj,, = k||+g'±[g^,-(k||+g')T/'e„ (6) 

where is the unit vector along the z axis, k|| is the reduced wave vector lying in the surface 
Brillouin zone (SBZ) of the given lattice, and g' is one of the reciprocal lattice vectors in 
Eq.®. 

In order to avoid the cumbersome formulas in the representation of the displacement 
field, we use the same complete set of spherical-wave solutions as recommended by Sainidou 
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et aliS, 



JLZ^(r) = ^V[ji{qir)Yrii)], Humi^) = ^V[/i+(gzr)rr 

JMimir) = ji{qtr)Ximir), HM/^(r) = /i+(gtr)X,„(f), (7) 

Jmmir) = X ji{qtr)Xi^{f), Hmm{r) = x /i+(gir)Xi„(r), 

here and are the spherical Bessel and Hankel functions, 1^"* and X;^ are the usual and 
vector spherical harmonics, respectively. 

The scattered waves by all spheres in the plane can be separated into two parts. One 
is the contributions from spheres of type A, u^'^(r), and the other is those from spheres of 
type B, u^'^(r'), here r = r' + ab. Let a^^^{P = L, M, N) be the scatter wave coefficients of 
sphere A at the site -R„. According to the Bloch theorem, we have 

u^/(r)= J2 <LH„„(r„)=5^a+f5^exp(zk||-R„)H„„(r-R„). (8) 

n,Plm Plm R„ 

where a'^{P = L, M, A^) are the scatter wave coefficients of sphere A at the origin. In the 
domain close to surface of sphere of type B, it is always satisfied that |r'| < |R„ — ab\ since 
all spheres are not overlapped. So that the u^/(r) can be expanded into spherical waves 
about the sphere of type B, which is located at ab, as follows 

= exp(2K^',,, • ab) ^ dZ^Pimir'). (9) 

Plm 

Likewise, we have 



uf/(r') = exp(zKf,, • 2) b^f^UpUO 

n.Plm 

= exp(iK^'^, ■ ab) ^Hm^ exp(ik|| ■ R„)Hpzm(r' - R„) 



(10) 



Plm R„ 

it can also be expanded as follows 



<'{r') = J2^Zjpim{r). (11) 



Plm 

In addition, the external incident wave is 



u,„(r) = J2 CJp/™(r) = exp(^Kj,, ■ d>) J] aZ^Pimir'). (12) 



Plm Plm 
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Moreover, if we remove the term corresponding to sphere located at = in Eq. ([8]), 
the terms left are just the scattering waves come from all spheres of type A except for itself. 
We denote it by u^'^'(r), and expand it into spherical waves about the origin 

^f:'ir) = J2aZ3pUr). (13) 

Plm 

In the same way, after removing the term corresponding to R„ = in Eq. (ITO!) . we have 

uf:'{r') = exp(zKj,, ■ ^) ^ b^JpUr')- (W) 

Plm 

It means the scattering waves come from all spheres of type B except for the one located at 
ab. It can be shown from direct calculation (see Appendix A for the details) that 

'P _ qPP' +P' ,'P _ 'Sp qPP' r+P' 

"Zm ~ ^Hm;l'm'"'l'm' ^ "im ~ 2-^ ^ Hm;l' m' " I' m' ^ 

P'l'm' P'l'm' Mq\ 

'P _ ST -PP' h+P' rl'P - \^ n+P' 

Hm — 2-^ ^lm;l'm'"l'm' ^ ^Im ~ 2-^ ^lm;l'm'"'l'm' ■ 

P'l'm' P'l'm' 

According to MST, the wave incident on the sphere of type A located at the origin consists 
of three parts, the first party is the externally incident wave, the second part is the sum 
of all the waves scattered by spheres of type A except for itself, and the third part is the 
sum of all the waves scattered by spheres of type B. The coefficients of the incident wave 
and scattering wave is related to each other by the Mie matrix, so we have the following 
relationship 

^tm ~ '^a^,lm;l'm'i'^Um' + '^l^m' + '^l^m') ■ (16) 

P'l'm' 

Likewise, the wave incident on the sphere B located at ab also consists of three parts, the 
externally incident wave, the sum of all the waves scattered by spheres of type B except for 
itself, and the sum of all the waves scattered by spheres of type A, so we have 

h+P- V tpp' fn^P' -^h'P' -^H'p'\ (^7) 

"im — 2-^ ''-b;lm;l'm' ["'I'm' + '^Vm' + "Z'm' ) ■ ) 
P'l'm' 
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Combining Eqs. (fTSi) . (fT6|) and (fT7|) . we obtain 



E 

P'l'm' 



A A A mPP" qP"P' 

OpP'Oll'Omm' ''-a;lm;l"m"^h"m";l'm' 



P"l"m" 



P'l'm' 



Erppp" '::P"P' 
^ a;lm;l"m"^l"m";l'm' 



P"l"m" 



,+P' _ yPP' OP' 

"«'m' ~ -'- a;lm,;l'rn'"'l'm'i 



P'l'm' 



P'l'm' 



E 



bdmd"m"^l"m"d'm' 



(18) 



+ E 



P"l"m" 

A A A mPP" qP"P' 

Opp'Ou'Omm' b;lm;l" m" ^ h" m" -,1' m' 



P'l'm' 



P"l"m" 



,+P' _ rj^PP' OP' 

"i'm' ~ -'-b;lm;l'm'^l'm'- 



P'l'm' 



When the coefficients of the incident wave are given, equations f|T8|) determine the 
scattering coefficients a^^ and 6;"^ of spheres of type A and B. Consequently, the waves 
scattered from the plane of spheres of type A and B are also obtained in terms of Eqs. ([8]) 
and (fTOj) . respectively. 



We write the coefficients in the form of 



Jg'j,' 



(19) 



based on Eq. (j5]), where coefficients A^^ are given in Psarobas et al'' (Eqs. (3.4), (3.7) and 
(3.8)). Because of the hnearity of Eqs. f|T8l) . the coefficients a^^ and h^^ can also be written 
as follows 

'-^''-EAt';.'«v)Mg'.', (20) 



"Im 



Km — ^h>^.i'(^e:'u' )[ 

i' 

so that the Eqs. (fT8|) are reduced to 



mls'i'i 



(21) 



P'i'm 



+ E 
E 



PTm' 



P"l"m" 

nPP" 'z:P"P' 
- a\lm\l"m"' — 'l"m" d'm' 



P"l"m" 



"y^, ^PP'^ll'^mm' T^,lm;l"m"^U'm";l'm' 
1 

E 

E 

p 

+ E 

P'l'm 

ErpPP' .OP' /IV-S' N 

''- b;lm;l'm' ^I'm' ■,i'\'^S'u' ) ■ 



P'l'm 



nPP" ^P"P' 
''bdmd"m"^l"m"d'm' 



P"l"m" 

A A A \^ T^PP" oP"P' 

Opp'OwOmm' b;lm;l"m"^ ^l"m" ■,l'm' 

P"l"m" 



^Vm';i'\'^?.'y') 

U+P' \ _ rj^PP' .OP' /Ty-S' N 

-°/'m';j'l-'^g'!^'J ~ 2-^ -'-a;lm;l'm'^l'm';i'\'^g'u')^ 
P'l'm' 

^^:''.'«'.o (22) 

^l^m'-A'i^s'u') 



P'l'm' 
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Now, the scattered waves given by Eqs. ([8]) and (fTOj) are expressed as sum of plane waves 
as follows (see Appendix B) 

= 5^[M^Jg,exp(2K^,-r)ei, 

^' , _ (23) 
uf/(r') = exp(2K^',^, ■ ab) ^[<f ]^ exp(2K|^ • r')e, = ^ [wfj ^ exp (iK^^ ■ r)ei. 

Here u'^ is obtained in coordinates r', it is obviously that its expression is similar to m^. 
The total scattered wave then becomes 

u:,(r) = J2i[ui]y + [nfj^) exp(^K|, ■ r)e„ (24) 

gi 

where [tt^jgj and [uf,]gj are given by 

= EE^r-;^(Kg.)^£.(Kj.')M^v, (25) 

i' Plm 

[«fc]g. = E E AL,(K^j5+^,(Kf,,)Mg'.' exp (^(Kj,, - K^J ■ 2). (26) 

i' Plm 

The derivation is given in Appendix B. 

With these results, we get this 2D "complex" crystal's M-matrix elements 

= Sss'S^^'Sw + E ^LAKuMtXAK^') + exp (z«,., - K^J ■ 2)]. 

Plm 

(27) 

Using the algorithm recommended by Psarobas et al^^^^, as soon as the M-matrix is given, 
both the frequency band structure of an infinite crystal and acoustic properties of a slab 
of this crystal can be calculated in the same way. Here we only write down the Q-matrix 
elements 

Q^^,^^' = ^^P [^(K+, • + ■ d,)] , 

Qgi;g'.' = ^gtg'.' H^i^ ■ - Kg,,, ■ d.)] , ^^^^ 
Q^^^,^, = exp [-z(Kg, ■ d, - ■ d,)^ 

Qilg'i' = ^^gi^'i' exp [-^(Kg, ■ d; + Kg,,, • d^ 



III. NUMERICAL RESULTS AND DISCUSSION 

We demonstrate the applicability of our method by applying it to two specific examples. 
Firstly, we illustrate the 2D periodic structure of NaCl-type phononic crystal in figure 1, 
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here the nearest distance between the same spheres is defined as lattice constant oq, so 
that ai = ao(l,0,0) and a2 = ao(0,l,0) in Eqs. ([1]) and ([2]). Then we choose d; = d,. = 
ao(0.25, 0.25, \/2/A) in Eq. ([28])^°, and obtain a slab of NaCl-type phononic crystal of finite 
thickness associated with the (001) surface. We take ao = 1cm in the calculation. 



o — o — o — o 





FIG. 1: The 2D structure of each layer of (001) plane is given in the left graph, the lattice constant 
oq is the nearest distance between the same spheres, the distance between two neighbor layers is 
ao/V2 for NaCl-type phononic crystal (the right graph is the section perpendicular to the layer 
plane) . 



In figure 2, we show the transmittance curves of a slab of 8 layers parallel to the (001) 
surface of NaCl-type phononic crystal at normal incident for longitudinal wave (a) and 
transverse wave (b). The materials are: sphere of type A is lead sphere of radius 5* = 0.25cm, 
sphere of type B is lead sphere of radius S = 0.10cm, the host matrix is epoxy. The 
relevant parameters are, for lead: p = 11.6 x lO^kg/m^, ci = 2490m/s, Ct = 1130m/s, 
and for epoxy: p = 1.18 x lO^kg/m^, q = 2540m/s, Ct = 1160m/s. Figure 3 shows the 
corresponding complex band structure of the infinite NaCl-type phononic crystal associated 
with the (001) surface. AH the units are following the suggestion given by Sainidou et aU^, 
where cq = 2540m/s. 

Comparing the results in figure 2 and figure 3, we see that there is a frequency band gap 
approximate between 1.6 to 2.25 for the longitudinal wave, and also several band gaps for 
the transverse wave. It should be noted that the width of the gaps and the positions have 
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FIG. 2: The transmittance curves of a slab of 
8 layers parallel to (001) surface of NaCl-type 
phononic crystal at normal incident for longi- 
tudinal wave (a the upper graph) and trans- 
verse wave (b the lower graph). 



FIG. 3: The phononic band structure at the 
center of the SBZ of a (001) surface of NaCl- 
type phononic crystal corresponding to Fig. 2. 
The filled circles represent the longitudinal 
mode, the open triangles represent the trans- 
verse modes and the open circles represent the 
deaf band, respectively. 



all changed from those of the traditional "simple" phononic crystal, so that a lot of work 
has to be done for searching some new properties of this kind of phononic crystals. 

In figure 4 we show the reflectance curve of a slab of layers parallel to the same surface at 
normal incident for longitudinal wave. The materials are taken as follows: sphere of type A 
is lead sphere of radius S = 0.20cm coated with a 0.5mm layer of silicone rubber, sphere of 
type B is silica sphere of radius S = 0.10cm coated with a 1.5mm layer of silicone rubber, 
the host matrix is epoxy. The relevant parameters are, for lead: p = 11.6 x lO'^kg/m^, 
Ci = 2490m/s, Ct = 1130m/s, for silicone rubber: p = 1.3 x lO^kg/m^, ci = 22.87m/s, 
Ct = 5.55m/s, for silica: p = 2.2 x lO^kg/m^, ci = 5970m/s, Ct = 3760m/s, and for epoxy: 
p = 1.18 X lO^kg/m^, ci = 2540m/s, Ct = 1160m/s. We take such materials that the spheres 
of type A and B have different local resonance frequencies, which can be estimated in the 
approach of Liu et al^. 

When the NaCl-type phononic crystal is consisted of two kinds of local resonance materi- 
als, our results show that both of local the resonance frequency (see the peaks with amplitude 
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close to 1 that represents total reflectance) are still keep unchanged, which means that this 
kind of hybrid does not change the property of intrinsic local resonance. In addition, as the 
number of the layers increased, the band gap is broadened indicate by the broader peak. 
We also see that there is a small wiggle near 0.11. In fact, when we scan the frequency 
in more refined steps, we found that there is a very sharp resonance located at there. The 
corresponding frequency band structure is shown in figure 5. 




FIG. 4: The reflectance curves of a slab of dif- 
ferent layers parallel to (001) surface of NaCl- 
type phononic crystal at normal incident for 
longitudinal wave for local resonance materi- 
als. The solid, dash, dot and dot-dash line 
are for 8, 4, 2 and single layer, respectively. 



FIG. 5: The phononic band structure at the 
center of the SBZ of a (001) surface of NaCl- 
type phononic crystal corresponding to Fig.4. 



The examples here are concerning the normal incident wave for simplicity, one can also 
take any incidence angles at will. Our results obtained with an angular momentum cutoff 
Imax — 4 and 13 g vectors, the relative error of the results are about 10~^. To obtain the error 
estimates, we have also performed the calculation with larger cutoffs of angular momentum 
and compared the results with different cutoffs. In addition, we have also compared our 
results of identical spheres, which becomes a "simple" crystal in fact, to the results of the 
known programm of "simple" phononic crystal for confirming our code. 
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IV. CONCLUSION 



We have shown that, for a system of two kinds of non-overlapping elastic spheres with 
different mass density, Lame coefficients and radius forming a NaCl-type structure, one can 
calculate the phononic band structure of the infinite crystal accurately and efficiently, and 
one can also obtain the transmission, reflection, and absorption coefficients of elastic waves 
incident onto a slab of the materials of finite thickness using the method and formalism of 
the present paper. 

The programm used in this work is a direct extension of the programm package by 
Sainidou et al^-, we thank the authors of the package for providing this useful package. 
The work is supported by the National Natural Science Foundation of China under grant 
No.90103035, No.l0334020 and No.l0174041. 

APPENDIX A 

In this appendix, we give the detail derivation of equation (ITSl) . Since Eq. ([8]) and Eq. 
dn]) are exactly the same, we have 

^ (il^Jpim(r') = exp(-i(k|| + g') ■ ab) ^ ^ exp(zk|| ■ R„)Hp,„(r' + ab - R„), (Al) 

Plm Plm R„ 

using the fact that exp(ig' ■ R„) = 1, it becomes 

5^0«n^(r') = ^a+f 5^ exp(z(k|| + g') ■ (R„ - ^))H„^(r' - (R„ - d>)). (A2) 

Plm Plm R„ 

The vector spherical waves originated at ab - 
the following expression (for example, see Eq 
matrix is just the matrix A here.), 

H„^(r' - (R„ - ^)) = AS„.'(-(Rn - ^))Jp'r,n'(r'), (A3) 

P'Vm' 

SO that 

^Z[vm' = $^exp[^(k|| + g') ■ (R„ - d))] A^^;,„,(-(R„ - d))). (A4) 



- R„ and those originated at ab are related by 
. (B9) in Appedix B of Sainidou et al^^, whose 
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Moreover, in the case of "simple" crystal, Sanididou et aU^. remarked that the evaluation 
of the matrices ^[rn 'i'm' becomes the evaluation of a well-known matrix Z in the theory of 
low-energy electron diffractionr-^. Here the matrices "^fXi'm' ^^"^ ^^^^ evaluated in the same 
way except that the definition of matrix Z is modified as follows 

= 5Zexp[z(k|| + g') ■ (R„ - 2)] G;^;/w(-(Rn - qt), (A5) 

where 



Gimr,n.{-T-qt) = Y,Eim{l'm'-l"m")h^{qtT)Yi"'\-f), (A6) 

I'm' 




In order to evaluate matrix Z in eq. (]A5|) . we introduce matrix K whose elements are 

^L"''(g*) = ''M^^H + g') ■ Rn] G',„;r,n'(-H.n; (A8) 

which can be numerically calculated by using Kambe's method^^^^. Based on the following 
relationship^^ 

h+{qt\Rn - ^1) + ^b)= iv_m'r'-m"'(ah- qt)h+,{qtRn)Yi;;r"'i-Rn), (A9) 

l"'m"' 

where 

^im;i"m"i-r;qt) = YEinril'm';l"m")jiiqtT)Yi^'{-f), (AlO) 

I'm' 

and some straightforward but trivial algebra calculation, we finally get 

Ziy'iqt) = eM~^ih+ S') ■^b]J2{Eimil'm';l"m") 

^'^ _^ (All) 

^[hUq,ah)Yr' {ah) + YKi^;{qt)jL{qtah)Y^^ {ah)]}. 

LM 

After that, we obtain the matrix Z through matrix K, so are matrix S. Similarly we get 
another matrix S in Eq. lfTSl) . In addition, it should be remarked that those two matrix S 
are different in general, except for a few cases, such as NaCl-type. For the two matrices Vl 
in Eq. (fT5l) . we have 

^^S'n^' = E ^^P(^k|| ■ R0AL^;,„,(-R0. (A12) 
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These matrices also appear in the case of "simple" crystal, its evaluation is already 
available^. 

However, there are some tricks for several specific examples to calculate the S matrix in 
eq. flA4|) . Here we propose the NaCl-type for instance. Supposing the spheres of type A is 
identical with those of type B, the crystal becomes a "simple" crystal and forms a new 2D 
periodic structure, its lattice is denoted by now. It is quite easy to obtain that 



' — 'Imd'm' 



J2 exp[^(k|i + g') ■ K]A[^',^,{~K) - ^^P[^(k|l + g') • Rn]ASw(-Rn). 

(A13) 

Based on the computer programm of "simple" crystal, it is obviously that only little modi- 
fication is needed for dealing with the NaCl-type crystal. This is actually our original idea, 
there are also other types can use this idea, for example, when ab = ai/2 or ab = a2/2. 

The specific form of the above matrix elements are given by Sainidou et al^°. It should 
be remarked that S have the same property of Q if the spheres of type A and type B locate 
in the same plane, it says 



^Im^Vm' ~ ^Irn^l'm' ~ ^hn;l'm' ~ ^5 UulcSS I + 171 + I' + m' : CVeU, 

Vm' = vm' = 0, unless l + m + l' + m': odd. 



(A14) 



APPENDIX B 

Firstly, we write down an identity^, 

exp(^kij ■ K^)hi{q^r^)Yr{rn) = ^/'"(K^J exp(zK± ■ r), (Bl) 

where Aq denotes the area of the unit cell of the lattice given by Eq. ([1]), the plus (minus) 
sign on Kg,^ must be used for z > {z < 0). 

Using the identity above, we write Eq. (jS]) as follows, 

= E E <^imi^'su) exp(^K^. ■ r), (B2) 

g Plm 
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where 



2TT(-iY 

^ qtAoK+,^l{l + l) I y ^t' 

+«r cos ^e-*F-+i (K^ J] e2 + * [ar'-e^*!^'"-' (K| J - 

qtAoK^t,y/l{l + 1) 

(B3) 

where a[" = a/(/ — + m + l)/2, and {6, (f)) denoting the angular variables of K^^. Sub- 
stituting from Eq. (!20!) into Eq. (\B2h . we obtain 

u^/(r) = J2[ui]l,expitKl, ■ r)e„ (B4) 

gi 

where 

[«^c]g. = EE^-;^(I^g-)At^.'(K^.OM^, (B5) 

i' Plm 
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